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Abstract 

We present magnetization measurements in single crystals of the tetragonal YN12B2C com- 
pound, which exhibit the phenomenon of peak effect as well as the second magnetization peak 
anomaly for H > 0.5 T (H || c). At the lower field (50 mT < H < 200 mT), we have observed 
the presence of flux jumps, which seem to relate to a structural change in the local symmetry of 
the flux line lattice (a first order re-orientation transition across a local field in some parts of the 
sample, in the range of 100 mT to 150 mT). These flux jumps are also observed in a single crystal 
of LUN12B2C for H || c in the field region from 2 mT to 25 mT, which are compatible with the 
occurrence of a re-orientation transition at a lower field in a cleaner crystal of this compound, 
as compared to those of Y N12B2C . Vortex phase diagrams drawn for H || c in LUN12B2C and 
YN12B2C show that the ordered elastic glass phase spans a larger part of (H, T) space in the 
former as compared to latter, thereby, reaffirming the difference in the relative purity of the two 
samples. 
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I. INTRODUCTION 



A well documented and researched issue in the context of vortex phase diagrams 1 ^^ of 
both low Tj*-&Z& as well as high T c superconductors 4 -"^ 1 ^ 1 ^ 1 ^ is the phenomenon of peak 
effect (PE) 2 , which is an anomalous increase in the critical current density (J c ) prior to 
reaching the normal-superconductor phase boundary. The PE is widely considered to mark 
a first order transition from a collective pinned ordered vortex solid to an individually pinned 
disordered solid 2 -* 1 ^ 4 - 11 ^. Another anomalous feature seen deep in the mixed state of weakly 
pinned samples of low as well as high T c superconductors is the second magnetization peak 
(SMP) anomaly^^ 1 ^. The SMP anomaly is often related to the pinning induced transition 
from dislocation free elastic glass (Bragg glass 3 ) to the dislocation mediated multi-domain 
vortex glass state. In such a circumstance, the onset field of the SMP anomaly is not ex- 
pected to vary with temperature 3 - 1 ^. 

Abrikosov2£ had predicted the flux lines to be arranged in a regular array, he found the pe- 
riodic array to be a square, but the difference in energy between the square and a triangular 
array is only 2 %. Numerous reports of change in the symmetry of the VL from a rhom- 
bohedral towards a square symmetry 2 ^ observed by small angle neutron scattering (SANS) 
measurements and Bitter decoration studies^ 1 ^ 3 - 1 ^ 4 ^ 1 ^ 1 ^ 1 ^ in the quaternary borocarbide 
compounds 29,30 , have resurged the interest in exploring and understanding the underlying 
mechanism governing this process in a wide variety of superconductors, e.g., PbTl 2 -, V^Sa 3 ^, 
La^ x Sr x Cu0 4 32 ^, Nb 34 , YBa 2 Cu 3 0^ etc. 

The borocarbide superconductors which have a tetragonal crystal structure with c/a ~ 3, 
are convenient test beds to study the interesting phenomenon of the change in the flux line 
lattice symmetry and the effects dependent on it, as high quality single crystals of large 
enough sizes can be grown by different procedures^. In the non-magnetic members of the 
borocarbide series, viz., YNi 2 B 2 C (Y1221) and LuNi 2 B 2 C (Lul221), it is observed that for 
H || [001] at low fields (~ 0.1 T), the VL symmetry is a distorted triangle (lower field rhom- 
bohedral Rl) with an apex angle j3\ < 60°. With the increase in field, the VL undergoes 
a sudden (first order) transformation via a 45° re-orientation to higher field rhombohedral 
Rh, with an apex angle j3 2 (> 60°). The Rh symmetry subsequently smoothly proceeds 
to a square symmetry via a continuous (second order) transition at a field H^^-. In the 
typical crystals of Y1221, the H x field (for H || [001]) lies in the range ~ 100-150 mT^ 4 ^, 
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whereas in the very clean crystals of Lul221, the same field lies in the lower field interval, 
20-50 mT for H || c 25,2 ^. Further studies^ 1 ^ in crystals of LUN12B2C have revealed that for 
H || [100], the rhombohedral Rh (apex angle fa > 60°) undergoes a sudden reorientation 
transition at H tr ~ 300 mT, such that the body diagonal of distorted rhombohedral locks 
up in [010] direction with apex angle ~ 82°. 

Kogan et al— predicted the occurence of the discontinous (first order like) Rl — > Rh tran- 
sition at ~ 20 mT for H 1 1 [001] in clean crystals of borocarbides. Subsequent studies by 
Gammel et alM in Co doped Lul221 crystals revealed that the Rh to square transition for 
H || [001] shifts from lower to higher fields with the increase in the Co doping, i.e., with the 
increase in disorder effects. However, the possibilities of relationship (s) between the struc- 
tural transition(s) in the VL and the spatial order-disorder transitions a la PE/SMP have 
not been described. It is of interest to know (i) how the symmetry transformations adjust 
to the pinning landscape, and (ii) whether the domain volume within which VL remains 
correlated depends on the underlying symmetry and the disorder effects on it. Magneti- 
zation studies by Silhanek et a/.— on a single crystal of Y1221 with an applied field H | 
c (i.e., [001]) revealed the presence of a kink in the pinning force density at a field value, 
which is close to the Hi value in this compound. To our knowledge, no other signature(s) of 
the VL symmetry transitions have been reported in the magnetization hysteresis measure- 
ments, while there have been several reports of the observation of the PE in the samples of 

Y1221 23 1 41 1 42 1 43 1 44 and Lul22 l™. 

We report here on the observation of the PE, the SMP anomaly and the flux jumps in the 
same isothermal magnetization hysteresis scan for H 1 1 c in the crystals of Y1221 and Lul221. 
The flux jumps interestingly occur in the field regime, where the VL symmetry transition is 
reported to occur across the respective Hi values for applied field oriented in the c-direction 
in each of the compounds. We believe that these flux jumps indeed have a correlation with 
the local symmetry change in the vortex lattice of borocarbide superconductors. To cor- 
roborate this possibility, we have traced in several ways, the minor magnetization curves by 
changing the initial thermomagnetic history of the sample of Y1221. The measurement of 
the quadrupolar signal, which purports to preferentially fingerprint the inhomogeneity in the 
magnetization across the sample, also, registers the change in the symmetry of the VL. The 
loci of the threshold fields at which the SMP anomaly and the PE commence can lead to 
the demarcation of boundaries across which changes in the spatial correlation of VL occur. 
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We present construction of the vortex phase diagrams for H 1 1 c in the crystals of Y1221 and 
Lul221. The parametric region in (H, T) space over which elastic glass state exists, seems 
to be influenced by the purity of the crystal. In Y1221 crystal, elastic glass state spans over 
a smaller region, indicating stronger pinning effect in this sample as compared to that in 
the crystal of Lul221. 

II. EXPERIMENTAL 

Magnetization measurements have been performed on two single crystals of Y1221, la- 
beled as A and B, and a crystal of Lul221. The single crystals, A and B, of Y1221 were 
grown by the travelling solvent floating zone method^, while the single crystal of Lul221 
was grown by the flux method, using N12B as flux 22 . The crystal A of Y1221 and that 
of Lul221 are (thin) platelets in shape, with the c-axis perpendicular to the plane of the 
platelet. The crystal B of Y1221 is, however, a parallelopiped (of size ~ 3 mm (1) x 0.7 mm 
(b) x 0.67 mm (t)) in shape, with the a-axis along the largest dimension. Both the crystals 
of Y1221 have T c (0) ~ 15.1 K, whereas the crystal of Lul221 has a transition temperature, 
T c (0) ~ 16.1 K. The DC magnetization measurements were performed using (i) a 12 Tesla 
Vibrating Sample Magnetometer (VSM) (Oxford Instruments, U.K.) and (ii) a 7.5 Tesla 
SQUID magnetometer (Model MPMS7, Quantum Design Inc., U.S.A.). 

III. RESULTS AND DISCUSSIONS 

A. Magnetization hysteresis measurements in YN12B2C 
1. Peak effect and second magnetization peak anomaly 

The PE in a weakly pinned crystal A of Y1221 is clearly evident (see Fig. 1(a)) by the 
characteristic 'bubble' in a five quadrant M-H loop recorded at 2.1 K using a VSM with a 
scan rate of 0.35 T/min and with an applied field H || c. Also, evident in Fig. 1 (a) is the 
presence of flux jumps at fields less than 200 mT, far below the peak effect region. The field 
values at which the flux jumps occur, vary somewhat from scan to scan, recorded at the 
same ramp rate of the swept magnetic field. This phenomenon will be described in detail in 
section B. 
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FIG. 1: A portion of the five quadrant M-H loop at 2.1 K in Y1221 (# A) with H || c. The presence 
of a peak effect 'bubble' and flux jumps (four in the third quadrant and three in the fifth quadrant) 
can be noted in the panel (a). Panel (b) shows a blow up of the M-H loop highlighting the second 
magnetization peak anomaly. The locations of the H™, H p and H C 2 are marked in this panel. 



Fig. 1(b) shows the M-H plot on an expanded scale to emphasize the presence of an anoma- 
lous feature designated as the second magnetization peak (SMP) anomaly around a field of 
1 T. The large hysteresis in the magnetization at fields less than 0.5 T in Fig. 1(a) gives 
way to a much smaller irreversibility at higher fields (upto close to the onset field H° n of the 
PE), indicating that the pinning in the sample is weak and an ordered elastic glass phase 
(i.e., a Bragg Glass (BG)^ like state) gets established well before the field regime of the 
PE. From the hysteresis plot of Fig. 1(a), it may seem that an ordered BG phase extends 
all the way up to onset field of the PE, but a closer look at the data in Fig. 1(b) (which 
reveals the presence of a SMP anomaly), suggests that the BG phase could terminate at 
the onset field of the SMP anomaly above which the multi-domain vortex glass (VG) phase 
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ensues^. The hysteresis width, AM(H) (oc J C (H)), starts to decrease with field once again 
above the peak field of the SMP anomaly and continues upto the onset field of PE, thereby 
implying an improvement in the spatial order in this interval. A complete amorphization of 
the multi-domain elastic VL commences only at H° n . 




Ho H ( T ) 



FIG. 2: Log-linear plots of normalized critical current density, J C (R)/J C (H P ) ( = AM(H) / AM(H p )) 
vs. field (/UoH) in Y1221 A) at two temperatures (2.1 K and 6.5 K) for H || c. The locations 
of the maximum position of the SMP anomaly (H smp ) and the onset (H° n ) and maximum (H p ) 
positions of the PE are marked. 

Fig. 2 shows a (normalized) J c versus H plot corresponding to the data of Fig. 1(a). Also, 
shown for comparison, is a (normalized) J C (H) plot (H || c), at a higher temperature ( T = 
6.5 K ). The peak fields of the PE (H p ) and the SMP (H smp ) stand marked for both the 
temperatures. H p decreases considerably with an increase in temperature (from 5.4 T at 
2.1 K to 3.7 T at 6.5 K), but H smp remains near 1.1 T at both the temperatures. This leads 
to an inference that the H smp (T) line in Y1221 A) would exhibit a very weak temperature 
dependence, reminiscent of the behavior of the SMP anomaly in another weakly pinned low 
T c superconductor, viz., Ca^RhiSrii^ (T c (0) ~ 8.2 K)^ZJ£ 

Fig. 3 shows a log-log plot of the current density (normalized to its peak value across 
the PE, where the correlation volume of VL is expected to reach a minimum) with field 
(normalized to the upper critical field H^) in Y1221 (# A) with H || c and at 6.5 K. Also, 
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FIG. 3: Log-log plot of normalized critical current density AM/AM p (oc j c ) vs. normalized field 
(h=H/iT c2 ) in Y122I (# A) and LuI221 at T = 6.5 K for H || c. The peak effect as well as the 
second magnetization peak anomaly are evident in the data of both the crystals. While only one 
power law region is evident for the case of Lul221, two power law regions, viz., PL-I and PL-II can 
be seen for Y1221. 

shown for comparison in this figure, is a plot in Lul221, at (nearly) the same temperature 
with H || c. Both compounds exhibit SMP like anomaly distinct from the PE feature, 
but with one little difference. While for Y1221, the onset field of SMP anomaly lying 
deeper inside the mixed state is well separated from the PE feature and does not display 
temperature variation, for Lul221, the SMP anomaly lies at the edge of the PE and varies 
with temperature as the PE. 

The field regime where j c varies with field in a power law manner, i.e., j c oc h~ n , is often 
demarcatec?^ as the collective pinned elastic regime. In Lul221, the power law behavior 
can be seen (cf. Fig. 3) to extend upto h (= H/H C 2) ~ 0.6, after which the BG phase 
probably gets broken up into multi-domain VG phase, which starts to amorphize at the 
onset of PE (h ~ 0.85). On the other hand, for Y1221, one can mark out two regions for the 
power law behavior. The first power law region (PL I) extends only upto h ~ 0.18 and the 
second region (PL II) surfaces between h ~ 0.4 and 0.75, where the dislocations injected in 
the interval, h ~ 0.18 - 0.4, could partially heal. 
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The ratio of J c at the peak position of the PE to that at the onset of the PE is about five 
for Y1221, while the same ratio is about two for Lul221. This implies that the correlation 
volume in Y1221 shrinks to about 1/25 of its value at the onset of PE (J c oc l/y/Vc), while 
for Lul221, V c at H p has shrunk only to 1/4 of its value at H° n . This is plausible, since for 
Lul221, the process of reduction in V c starts at ff°" p and it continues till the arrival of H° n . 
On the other hand, for Y1221, the process of reduction in V c starts at the onset of PE. Prior 
to it, the VL in this case has a possibility to heal between H smp and H° n . In fact, the ratio 
of J c at H p to its value at H°^ p (h ~ 0.6) is also about five in Lul221. The premise that 
the SMP anomaly is disorder induced could imply that its onset at h ~ 0.18 in the given 
Y1221 crystal A) as compared to the onset of SMP at h ~ 0.6 in Lul221 crystal, reveals 
the relative purity (levels of effective disorder) of the samples of these two compounds. 

2. Pinning force density 

Fig. 4 shows a plot of the pinning force density (oc J c x H) with field on a linear-log 
scale in Y1221; the data correspond to the magnetization hysteresis measurements shown in 
Fig. 1(a). Three peaks, marked as (F p ) max , SMP and PE, can be clearly distinguished. An 
increase in the pinning force upto (F p ) max is a representation of the increase in the rigidity 
of the vortex lattice as a result of the interactions between the flux lines. A peak at (F p ) max 
probably implies that the rigidity of the vortex lattice has attained a limiting value dictated 
by the interaction effects. The second hump corresponds to the SMP anomaly and the third 
peak corresponds to the quintessential PEA&. Apart from the three maxima, one can also 
mark out the positions of three flux jumps in this figure. These jumps are observed to lie 
close to, though a little lower, the field value corresponding to the (F p ) max . This suggests 
that flux jumps are observed in the field regime, where the lattice rigidity is effective and 
VL is well formed. 

B. Flux jumps in yA^-E^C 

1. Zero field cooled measurements in crystals of Y1221 

Panel (a) of Fig. 5 shows a five quadrant M-H loop recorded at 2.1 K in a VSM at a sweep 
rate of 0.25 mT/sec in Y1221 (# A) with H || c. One can clearly discern the presence of the 

8 



°' 3 1 YNi 2 B 2 C(#A) 

H || C I 

_ T = 2.1 K 

w 

c 0.2 - 

3 

S3 
k. 

CO 

F Flux 




H H(T) 



FIG. 4: Plot of the pinning force density (J c xH) vs. field in Y1221 (# A) with H || c. The 
regimes of (F p ) ma:r , SMP and PE have been indicated, and the positions of the flux jumps have 
been marked by arrows. 

multiple flux jumps in selective quadrants. The following observations are noteworthy: 

(i) Flux jumps are absent in the ZFC run (0 — > H max ), in the second quadrant (H max — > 
) and in the fourth quadrant {-H max — > 0). They occur only in third (0 — > —H max ) and fifth 
(0 to H max , subsequent to the initial ZFC run) quadrant. 

(ii) Two flux jumps occur at about -55 niT and -120 mT, in third quadrant and three flux 
jumps happen at 70 mT, 110 mT and 150 mT in fifth quadrant. 

On cycling the field repeatedly between ± 300 mT several times, it was noted that the 
number of flux jumps in third and fifth quadrants stabilized to three (see panel (b) in 
Fig. 5), however, the precise field values at which the jumps happen were found to vary each 
time. In order to overcome an apprehension that these flux jumps could be an artefact of the 
rapid ramping of the magnetic field and/or the procedure of magnetization measurement in a 
VSM, M-H loop was also recorded using a SQUID magnetometer, where the superconducting 
magnet is kept in the persistent mode while ascertaining the magnetization value of the 
sample. Qualitatively, the same behaviour (data not shown here) as depicted in Fig. 5 was 
noted. The number of flux jumps in one of the quadrants (namely, third) were, however, 
observed to increase from three to four. We are inclined to surmise that the number of 
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FIG. 5: Panel (a) shows a five quadrant M-H loop in Y1221 (# A) at T = 2.1 K with H || c using a 
VSM with a scan rate of 0.25 mT/sec. Presence of two flux jumps in the third quadrant and three 
flux jumps in the fifth quadrant can be clearly noted. Panel (b) shows the M-H curve obtained 
after cycling the field several times between ±H max (= ± 300 mT). Note that the number of flux 
jumps in third and fifth quadrant is now same. 

jumps in the third/fifth quadrant could statistically vary between two to four. 

To establish the notion that the flux jumps do not depend on any specific physical char- 
acteristic of a given sample of Y1221, the M-H loops were recorded on another sample of 
Y1221 (crystal B, which is parallelopiped in shape) with field applied parallel to its longest 
physical dimension (H || a in this case, see Fig. 6). The presence of flux jumps in the third 
and fifth quadrants attests to the fact that demagnetization factor of the sample does not 
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influence the manner in which the flux jumps get observed. The observation of flux jumps 
for H 1 1 a orientation in Fig. 6, also, suggests that the mechanism responsible for these jumps 
probably does not depend on the orientation of the applied field w.r.t. the crystalline axis. 

It may be pertinent to state here that flux jumps in the third and fifth quadrants have been 
noted in the hysteresis loops pertaining to the intermediate state of some of the specimen 
of Type-I superconductors, viz., Pb 46 . It is believed that in such a circumstance, the flux 
jumps are caused by the escape of flux by the fusion of local macroscopic regions of positive 
and negative magnetization lying in juxtaposition. 

In order to comprehend the mechanism governing the flux jump process, field profiles relevant 
to Fig. 5 are drawn, as per a simplified prescription of the Bean's Critical State model 4 - 
for an infinite slab of thickness 2D, with surfaces at x = ± D, the applied field H being 
parallel to the surface. We realize that our experimental sample is either platelet shaped 
or is a parallelopiped, but it is hoped that the Bean's profile drawn in Fig. 7, for zero 
demagnetization limit and for J c (H) = constant for the sake of simplicity, would qualitatively 
turn out to be instructive even when J c (H) decreases with an increase in H. For brevity, 
profiles are shown in only one half of the sample. The other half is a mirror image of the 
drawn profiles. While drawing Bean's field profiles, B-x (local macroscopic field, B, versus 
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FIG. 7: Bean's profiles (only one half section (-D to 0) shown) relevant for the five quadrant M-H 
loop in Fig. 5 (a). 



distance, x, from the centre of the sample), we have assumed that the full penetration field, 
H*, is 200 mT. In the prescription of Bean's model, H* is the limiting field at which the 
applied field invades the entire sample, after initial zero field cooling the sample (see panel 
(a) of Fig. 7). The estimate of ~ 200 mT for H* in Y1221, with H 1 1 c, at 2.1 K is based on the 
magnetization data in Fig. 5 (a). It is the limiting field where the virgin ZFC magnetization 
curve would merge with the forward leg of the envelope^ loop in the absence of the flux 
jumps. When J c (H) decreases with field, such a limiting field can be taken to identify the 
nominal H* value. 

In Fig. 7, we have restricted ± H max to ± 300 mT in conformity with the M-H data shown in 
Fig. 5(a). Panels (c) to (e) of Fig. 7 are relevant for the M-H loop shown in panel (b) of Fig. 5 
as well. It may be noted that the field profiles in panels (a), (b) and (d) of Fig. 7 are such 
that the magnetic field inside the sample remains of the same sign (positive or negative) in 
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these cases. The field profiles in panels (c) and (e) of Fig. 7 allow for the possibility of both 
positive and negative field values inside the sample. The positive and negative (macroscopic) 
fields may be identified with the domains (or local regions) comprising vortices and anti- 
vortices. A narrow region encompassing zero field value would be free of any kind of vortices, 
and in its neighborhood, the domains of vortices and anti-vortices would lie in juxtaposition. 
The fact that the flux jumps are observed in the third and the fifth quadrants in conjunction 
with the Bean's profiles in Fig. 7 could further imply that co-existence of domains of vortices 
and anti- vortices is necessary for the jumps to occur. 

We would like to now surmise that mere juxtaposition of the regions of vortices and anti- 
vortices is not adequate to trigger flux jumps. We conjecture that sudden annihilation of 
vortices (flux jumps) gets triggered, when a change in the local symmetry of the vortex lattice 
from low field rhombohedral Rl to higher field rhombohedral Rh (or vice- versa) occurs over 
a portion of the sample (i.e., in some domains), while the regions of positive and negative 
fields lie in juxtaposition anywhere else in the sample. Let us call the domains comprising 
rhombohedral Rl and negative(positive) field as R£ , and, similarly, the domains with 
rhombohedral Rh and negative(positive) field as R H ■ The R^(R£) — > R H (R H ) (and 
vice-versa) transition(s) can be expected to happen in Y1221 crystal as the (macroscopic) 
field inside the sample crosses the region of 100-150 mT during the ramping of the external 
field. In addition, the Rh domains can, however, exist in (metastable) supercooled state 
during the field cooling process even in a field value somewhat less than 100 mT, as has been 
reported in some of the SANS experiments^ 1 ^. Such metastable domains could display the 
tendency to suddenly transform to the stable domains, during the subsequent field ramping 
cycle. 

On the basis of the above conjecture, let us now reexamine the observations in Fig. 5. During 
the first ramp up of the field after ZFC (0 mT to 300 mT), the disordered bundles of vortices 
will enter the sample from the corners and edges and attempt to settle down into regions of 
Rl and R H domains, with the former lying near the centre for H = 300 mT (see Fig. 7 (a)). 
As the applied field is reduced to mT (see Fig. 7 (b)), the R~l~ domains would be present 
near the sample edge and R H domains are likely to lie in the interior. No flux jumps happen 
in the first two quadrants, as the domains with anti- vortices are not present. In the third 
quadrant, as the field gets cycled to negative values, the first flux jump happens near about 
-50 mT. The field profiles in Fig. 7 (c) imply that in such a circumstance, the R~[ domains 
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would lie near the edge, while R\ and (supercooled/metastable) Rjj domains could exist 
in the interior of the sample. A change in the symmetry from Rjj — > in some domains 
could trigger an additional perturbation in the dynamically varying distribution of R~[, R^ 
and Rft domains, such that the vortices and anti- vortices lying in juxtaposition in R^/R^ 
domains could start annihilating each other leading to an avalanche resulting in a flux jump. 




-100 
-200 



100 
70 



(b) 




FIG. 8: Bean's profiles (only one half section (-D to 0) is shown) at an applied field of hqH = -100 
mT corresponding to a situation (a) just before the flux jump and (b) just after the flux jump, (c) 
Bean's profiles as the field ramps upto ^oH = -140 mT 

From Fig. 5, one can find that the magnetization values start to build up once again 
after the flux jump happens, and one can envisage the notion of an underlying envelope 
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M-H loop. As the applied field ramps up beyond -100 rriT, another (more) flux jump(s) 
happen (cf. Fig. 5 (b)). Such jumps could possibly be triggered by R~l — > transition in 
the domains, which are nucleating near the edge of the sample. When the applied field is 
in the range of -100 mT to -150 mT, the R\ domains would be present in the middle of 
the sample. Their presence would permit the annihilation of vortices and anti-vortices in 
R\ and R~l domains lying in juxtaposition in the interior. The recurrence of the process 
of R~l — > Rjj transition could trigger more than one avalanche, while the field is getting 
ramped down to -200 mT in the third quadrant. Similar arguments will explain the flux 
jumps in the fifth quadrant. 

We show in Fig. 8 a plausible sequence of Bean's profile, when the flux jump is 
located, say at H = 100 mT. As the flux jumps occurs, the Bean's profile of Fig. 8 (a) 
probably transforms to the profile in Fig. 8 (b), due to movement of anti- vortices from the 
exterior of the sample into the portion in which anti-vortices and vortices have annihilated 
to create a current free region at the right side of B = 0. The profile in Fig. 8 (b) implies 
the rearrangement of vortices in the interior of the sample, as shown by the solid line. It 
is apparent that profile in Fig. 8 (b) corresponds to a lower net magnetization value as 
compared to that for profile in Fig. 8 (a). As the field would further ramp away from - 
100 mT, the critical current would gradually get set up in the gradient free region. Fig 8 
(c) shows that at -140 mT, the Bean's profile has assumed a form, as if the flux jump had 
not occurred at -100 mT, thereby implying the return of the magnetization values to the 
underlying envelope hysteresis loop. 

2. Tracings of the minor hysteresis curves with different thermomagnetic histories 

To check the validity of our conjecture, we traced several minor hysteresis curves with 
different thermomagnetic histories, including the tracings of the complete hysteresis loops 
after having cooled the sample in different fields. Fig. 9 shows two representative minor 
hysteresis curves obtained at 2.1 K in Y1221 (# A) using the field ramp rate of 0.2 T/min 
in the VSM. In panel (a) of Fig. 9, after zero field cooling, the field is initially ramped upto 
+75 mT (filled squares), it is then reversed to -500 mT (open triangles), followed by ramping 
up again to +75 mT (filled circles). No flux jump is observed at any magnetic field in the 
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third quadrant, however, a flux jump occurs at about +32 mT in the fifth quadrant. 
Panel (b) of Fig. 9 shows that if the field is reversed to -500 mT from +174 mT, two flux 
jumps are observed at about -91 mT and -137 mT, respectively, in the third quadrant. It is 
pertinent to note that a flux jump in the range of -30 mT to -50 mT is not present in the 
third quadrant. On ramping up the field from -500 mT to +174 mT , one can witness three 
flux jumps at about 37 mT, 97 mT and 142 mT, respectively, in the fifth quadrant. 



YNi 2 B 2 C(#A) 




H H (mT) 

FIG. 9: In panel (a), the filled square data points trace the initial magnetization curve (0 to 75 mT) 
after zero field cooling (ZFC) at T = 2.1 K in Y1221 (# A) with H || c. The open triangle and 
filled circle data points in different panels correspond to the minor curves traced while ramping 
up the field to -500 mT from (a) 75 mT and (b) 174 mT respectively, and reversing the field from 
-500 mT to (a) 75 mT and (b) 174 mT, respectively. 

Fig. 10 shows the Bean's profiles corresponding to the data in Fig. 9. The profiles in 
Fig. 10 (a) imply that when the applied field reaches a value of -75 mT, the R£ and R~£ 
domains lie in juxtaposition in the outer portion of the sample and there are no vortices of 
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(a) Forward from to 75 mT 
and return to -500 mT 



(c) Return from 
174 mT to -500 mT 



174- 




FIG. 10: Panels (a) and (b) correspond to Bean's profiles (only one half section (-D to 0) shown) 
relevant for the minor curves of Fig. 9 (a). Panels (c) and (d) show the same for the minor curves 
of Fig. 9 (b). 

any kind in the middle region of the sample. As the applied field gets ramped from -75 mT 
to -150 mT, the specimen gets progressively more filled with R£ domains. As the ramping 
of the field proceeds further from -150 mT to -500 mT, the R~jj to transition could 
(repeatedly) occur near the edge of the sample, but there are no domains with vortices of 
the opposite kind like or Rjj like) in the interior. The conjecture stated above thus, 
precludes the occurrence of flux jump in the third quadrant in Fig. 9 (a). The profiles in 
Fig. 10 (b) further show that when the flux jump happens on approaching +37 mT field in 
Fig. 9 (a), the R£ and R~l domains would lie in juxtaposition near the edge of the sample 
and R]j domains would lie in the deep interior, near the central region. A transition from 
R~fj to R~l in the interior region probably provides the necessary trigger for annihilation of 
vortices across the boundaries of R~£ and R^ domains. 

An examination of the profiles in Fig. 10 (c) in conjunction with the flux jumps in Fig. 9 
(b) would imply that as the applied field reaches a value of -50 mT in the third quadrant, 
the R~l domain will lie near the edge, and only the R~l domains would fill the interior of the 
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sample. There is no possibility of Rl — > Rh transition anywhere in the sample between 
and -50 mT. The R~l — > R H transition probably happens near the edge of the sample as the 
applied field ramps towards -91 mT (see Fig. 9 (b)), at that stage, the interior of the sample 
contains R~l~ domains. As the field ramps further to -137 mT, another R^ to R H transition 
could occur near the edge, while the R~l domains are still left in the interior of the sample. 
This could explain the two flux jumps in the third quadrant in Fig. 9 (b). The profiles in 
Fig. 10 (d) further show that as the applied field approaches +37 mT in the fifth quadrant, 
the situation is similar to that in the fifth quadrant as depicted in Fig. 10 (b) and the first 
flux jump in Fig. 9 (b) gets triggered at such a field, by R H to R~[ transition in the interior. 
The subsequent two jumps in the fifth quadrant in Fig. 9 (b) are triggered by the possibility 
of recurring to R H transitions near the sample edge, while R~l domains continue to exist 
in the interior. 

Field cooled (FC) measurements, in principle, result in a near uniform field distribution 
across the sample. In a FC(H) state, the crystal will comprise Rl/Rh or a co-existence of 
Rl and Rh domains, depending on the field value. Tracings of minor hysteresis loops by 
changing the field in the third quadrant, from different Mpc(H) values revealed (data not 
being shown here) one/two/three jumps. Their occurence can also be rationalized in terms 
of above stated conjecture by drawing appropriate Bean's field profiles. 

3. Temperature dependence of the flux jumps 

It is well documented that the tendency and the magnitude of the flux jumps decreases 
with an increase in temperature^. The plots of the M-H loops in Y1221 (# A) at higher 
temperatures (T > 2.1 K) in Fig. 11 seem to conform to this notion. A comparison of the 
data in Fig. 11 with the corresponding data in Fig. 5 (a) reaffirms the premise that the 
extent of jump(s) decreases on progressive increase of temperature from 2.1 K to 8.45 K. 
A significant observation is that the flux jump evident at lower fields (50-75 mT range) in 
Fig. 5 (a) is not present at higher temperatures (i.e., T > 4 K as in Fig. 11). This trend 
could also find a rationalization in terms of the decrease of «/ c (H) (or full penetration field) 
with the increase in temperature, and the corresponding changes in the field distribution 
inside the sample. In Fig. 11 (b), the H* value at 6.45 K would be reckoned to be about 
150 mT. Bean's profile drawn for H* = 150 mT would then demonstrate the absence of flux 
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jumps at the lower field values. 
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FIG. 11: M-H loops at 4.45 K (a), 6.45 K (b) and 8.45 K (c) in Y1221 (# A) with H || c. Inset in 
panel (c) shows temperature variation in H reor and (see text for details). 

Our assertion that a flux jump is triggered by ^ ^ Rh transition and is facilitated 
by the juxtaposition of vortices and anti-vortices implies that the (highest) limiting field 
value at which a jump is observed (in a given fifth/third quadrant) has to be lower than the 
corresponding H* at a given temperature. Such a limit could also be taken as indicative of 
the field value at which — > Rh transition happens near the edge of the sample. Keeping 
this in view, we draw attention to the plot of such limiting values (designated as H reor and 
marked by arrows in panels (a) to (c) of Fig. 11) as a function of temperature in the inset 
of Fig. 11 (c). Multiple values of H reor at a given temperature represent the spread in these 
values during different M-H runs at the chosen temperatures. It is satisfying to see the 
similarity in }{ reor {T) line determined as above from our flux jump data and the i?i(T) line 
(representing R^ — > Rh transition) determined from SANS measurements in a crystal of 
Y1221 by Dewhurst et <M. 
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C. Flux jumps in LUN12B2C 
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FIG. 12: Expanded portions of the M-H loop obtained at 2.2 K in Lul221 with H || c. Panels (a) 
and (b) depict the M-H curves in the fifth (ramping the field from to H max = 200 mT) and the 
third (ramping the field from to -H max = -200 mT) quadrants, respectively. The locations of the 
flux jumps have been marked with arrows. 



The plots in Fig. 13 show parts of the M-H loop recorded using VSM at 2.2 K in Lul221, 
for H 1 1 c, with a scan rate of 0.25 mT/sec. The presence of the three flux jumps in the third 
and fifth quadrants in the field interval 2 to 25 mT can be noted. From Bitter decoration 
experiments performed in the field cooled state on single crystals of Lul221, the Rl — > Rh 
transition is reckoned to occur in the interval of 20 mT to 50 mT— . The observation of flux 
jumps in a single crystal of Lul221 at smaller field values as compared to those in Y1221 
is encouraging. The nominal H* at 2.1 K in the given crystal of Lul221, with H || c, is 
about 30 mT. Simplistic reasoning based on Bean's profiles for H* = 30 mT can rationalize 
the occurrence of the flux jumps at the observed field values. The local macroscopic field 
would envisage Rjj to R~l transition deep inside the sample and R£ to Rjj transition near 
the surface as the applied field ramps from zero to -30 mT in the third quadrant. 
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D. Measurement of quadrupolar signal (Q) in the crystal 'A' of YN12B2C using a 
VSM 



When the magnetization in a given sample is non-uniform, multipole moments other 
than the dipole, also, contribute 50 to a measured signal. Preferential measurement of the 
quadrupole moment (Q) is a very useful technique to gain information on the spatial in- 
homogeneity in the magnetization across a given sample 51 . Such a measurement has been 
performed using a VSM in Y1221 A), with H 1 1 c (see Fig. 14) to explore the fingerprints 
of the notion of the flux jumps observed in the M-H loops. The details of the measurement 
procedure can be found elsewher© 5 ^^. In brief, a sample is moved from the central region of 
the astatic pair of the coils of VSM to another location, where the signal due to the dipole 
moment is expected to cross from a positive to a negative value. At such a location, the 
measured signal preferentially captures the contribution from the quadrupolar moment of 
the sample. 

Panels (a) and (b) of Fig. 14 show Q (in arbitrary units) vs. H plots at 2.0 K and 2.9 
K, respectively as the field is sweeped from zero to H max in the so called fifth quadrant 
at a sweep rate of 0.01 T/minute. The Q vs. H data in Fig. 14 (a) appears to have 
a shape similar to the M vs. H data in Fig. 5, and this could imply that the residual 
contribution from the dipole moment of the sample is still dominating the measured signal 
at the new preferred location for the record of the quadrupolar contribution. There is, 
however, one notable observation in both the panels of Fig. 14, which purports to support 
the significant presence of the contribution from the quadrupolar moment of the sample 
at the new location. It may be mentioned that at the field values of the flux jumps, the 
measured signal gets out of the range of the plot with the Lock-in amplifier settings getting 
overboard momentarily and it returns to within the range only when the process of the flux 
jump is complete. A motivated search for such an occurrence (i.e., signal getting out of 
range) in the magnetization response at the usual central position of the astatic pair of the 
coils, did not yield an affirmative answer. This implies that the residual presence of the 
dipolar signal at the preferred position for the Q measurement, is not responsible for the 
signal getting out of range, when the flux jump happens. 

We are tempted to conjecture that a readjustment in the vortex matter during the flux 
jump process causes the contribution from the quadrupolar moment to undergo a peak 
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FIG. 13: The panels (a) and (b) depict the plots of quadrupolar signal Q (in arbitrary units) vs. 
H in Y1221 (# A) for H || c at T = 2 K and 2.9 K, respectively. 

like behavior. We believe that a study of anomalous variations in the quadrupolar signal, 
conveniently measurable in a VSM, has the potential to reveal the changes in the state of the 
vortex matter, which may not get fingerprinted as anomalous variations in the magnetization 
hysteresis response, like, the SMP anomaly and/or the peak effect phenomenon. 

E. Phase Diagrams in YN12B2C and LuNi 2 B 2 C for H || c 

Collating all the data together, we draw the phase diagram for Y1221 A) and Lul221 
for H || c. The diagram comprises H reor , H™ p , H° n and H c2 lines. The H™ p (T) lines 
have been drawn over a limited temperature interval over which data pertaining to it are 
presently available in these crystals. It can be noted that H^ p line for Y1221 is flat, while 
it decreases with temperature for Lul221. The former behaviour is reminiscent of the SMP 
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line in C a^Rh^Sni^^- and in optimally doped Bi2Sr2CaCu20s crystals^, while the latter 
is similar to the SMP line in 2H-NbSe2 crystal^. 
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FIG. 14: Vortex phase diagrams in YNi 2 B 2 C (# A) and LuNi 2 B 2 Cwith H || c. # reor , i?, 
i?p n and ff C 2 lines have been sketched and different phases have been named.. 
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The SMP anomaly is believed to mark a transition from a dislocation free elastic BG 
phase to a multi-domain VG phase-' 19 . The region below H°^ p and above H reor marks the 
BG phase. From Fig 15, it is also evident that the BG state spans over a larger (H,T) space 
in Lul221 as compared to that in Y1221, attesting to the fact that lesser residual disorder 
prevails in the former sample. 

IV. SUMMARY 

To summarize, we have expounded the details of magnetization hysteresis loops in single 
crystals of YNi 2 B 2 C for H || c. In the high field regime, a second magnetization peak 
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anomaly and the peak effect are observed, while at lower fields (H < 200 mT), the flux 
jumps are evident. The flux jumps have been noted in minor hysteresis curves recorded 
with differing thermomagnetic histories, in order to explore a probable reason for their 
occurrence. On the basis of Bean's profiles sketched for different thermomagnetic histories, 
we have conjectured that the flux jumps get triggered as a result of the structural transitions, 
Rl,h — ► Rh,Li possible in the vortex lattice in the field range 100-150 mT for H || c in 
YNi 2 B 2 C. Similar flux jumps have also been observed in a cleaner crystal (as compared 
to those of YNi 2 B 2 C) of LuNi 2 B 2 C at lower fields (H < 25 mT), giving support to an 
assertion made by us. Another evidence in support of our conjecture could come from the 
measurement of quadrupole signal (Q) in YNi 2 B 2 C using a VSM. Q undergoes a peak like 
feature at applied field values, corresponding to which the macroscopic field somewhere inside 
the sample is such that Rl,h —> Rh,l transition happens. The onset of a SMP anomaly at 
a lower field in the crystal of YNi 2 B 2 C as compared to that in LuNi 2 B 2 C could be termed 
as consistent with the occurrence of Rl — > Rh transition at a higher field in it. Based on 
all the measurements, vortex phase diagrams for Y1221 A) and Lul221 for H || c have 
been drawn, which depict various phases of vortex matter. Our results call for a study to 
map Bean's profiles by local micro-Hall bar arrays^ in all the five quadrants for different 
thermomagnetic histories, in samples which display flux jumps. 
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